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The Dyson-Schwinger equations (DSEs)B provide an approach well-suited to the cal- 
culation of pion observables. Since a chiral symmetry preserving truncation scheme 
exists,el they provide a framework in which the_dichotomous bound-state/Goldstone- 
mode character of the pion is easily captured.au Furthermore, because perturbation 
theory is recovered in theLweak coupling limit, they combine; e.g., a description of 
low-energy it-it scatteringEl with a calculation of the electromagnetic pion form fac- 
tor, Fn{q 2 ), that yieldsH the l/(j 2 -behaviour expected from perturbative analyses 
at large spacelike-g 2 and a calculated evolution to the p-meson pole in the timclikc 
regionE These and other contemporary applications are reviewed in Refs. H||. 

As an illustration, using the inhomogeneous Bethe-Salpeter equations for the 
axial-vector and pseudovector vertices; the dressed-quark DSE; and the fact that 
a nonperturbative Ward-Takahashi identity preserving truncation of the DSEs is 
possible, it was shown in Ref. j| that for flavour nonsinglet pseudoscalar mesons 

f H m 2 H = M C H r e H , (1) 

with M C H := tr flavour [M (c) {T H , (T fl )'}] , where M (f) = diag(m£, m% m$, . . .), with 
C the renormalisation point, and (•)* indicates matrix transpose, so that M\j is the 
sum of the constituents' current-quark masses. This model-independent identity is 
valid for all current-quark masses, irrespective of their magnitude, and therefore 
provides a single formula that unifies the light- and heavy-quark regimes. 
In Eq. (fy), fn is the leptonic decay constant, 



f H P, = Z 2 l — -tr 



T^^xh^P) , (2) 
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where Zi = ^(C, A) is the dressed-quark wave function renormalisation constant, 
with A the regularisation mass-scale. (The advantages of employing a translation- 
ally invariant regularisation scheme when Ward-Takahashi identities are involved 
should be obvious.) The r.h.s. is gauge-invariant and independent of £ and A. The 
other factor is 



i =Z,J^ itrfCr^XffteP)] , (3) 



where Z4 = Z±(£,A) is the dressed-quark mass renormalisation constant. Here 
the r.h.s. is gauge invariant and cutoff-independent, and Z^s ^-dependence ensures 
that the product on the r.h.s. of Eq. ([!]) is independent of the renormalisation 
point. Using these formulae it can be established that Eq. ([!]) reproduces the so= 
called Gell-Mann-Oakes-Renner relation in the limit of small current-quark massesB 
and also predicts tha±_heavy-meson masses increase linearly with the mass of their 
heaviest constituent However, the latter "heavy-quark limit" provides a poor 
approximation to the pseudoscalar meson mass trajectory for current-quark masses 
less than that of the 6-quark.liil 

Another more recent application is a calculation of the pion's valence-quark 
distribution,!^ of which a very accurate measurement is possible given a high- 
luminosity electron-proton collider.E^I Using an algebraic DSE-modelJi-3 employed 
successfully in studies of a wide range of hadronic observables; e.g., Refs. pO| , pT| , it 
is straightforward to calculate the "handbag" contributions to the virtual photon- 
pion forward Compton scattering amplitude. These are the only impulse approxi- 
mation diagrams that survive when calculating the pion's structure function in the 
deep inelastic Bjorken limit and yield the results in Fig. |l|. 

In this calculation* dressed-quarks with a valence-quark mass of M = 301 MeV 
carry 71 % of the pion's momentum at a resolving scale go = 0.54 GeV = 1/(0.37 fm). 
The remainder is carried by the gluons that effect the binding of the pion bound 
state. The second and third moments of the distribution are (x 2 ) qo = 0.18, (x 3 ) qo = 
0.10. To determine the resolving scale, go, we employed the 3-flavour (Aq^ 3 = 
0.242 GeV), leading-order, nonsinglet renormalisation group (evolution) equations 
to evolve the distribution in Fig. 0up to q — 2 GeV, and required agreement between 
the first and second moments of our evolved distribution and those calculated from 
the phenomenological fits in Ref. ]lq|. go = 0.54 GeV gives 
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0.10 ±0.01 
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0.27 ±0.01 


0.11 ±0.3 


0.048 ± 0.020 



(4) 



with the valence-quarks now carrying a momentum-fraction of 0.49. 

The evolved distribution is also shown in Fig. [l]. The evident accentuation via 
evolution of the convex-up behaviour of the distribution near x — 1 is characteristic 
of the renormalisation group equations, which populate the sea-quark distribution 

"Herein we have corrected a minor numerical error discovered in the calculations of Ref. . 
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Figure 1: Dashed line: calculated form of xu v (x; qo); Solid line: evolved dis- 
tribution, xu v (x;q = 2GeV); Dotted line: xu v (x;q = 4.05 GeV), evolved from 
xu v (x;q — 2 GeV) with Aqq^ = 0.204 GeV; and Dot-dashed line: phenomenolog- 
ical fit in Ref. Jl5| , which takes the form xu v (x) oc x a (l — x)@. The dataEl are 
obtained with an invariant /i + /i~-mass > 4.05 GeV and inferred from the differen- 
tial pion-nucleon Drell-Yan cross section using simple distribution parametrisations 
of the type just indicated, yielding a = 0.64 ± 0.03, (3 — 1.15 ± 0.02. This data was 
part of the set employed in the fit of Ref. |l^] . Fits to our calculation using the 
parametrisation in Eq. (^) are indistinguishable from our result on the scale of this 
figure. 

at small- a; at the expense of large- a; valence-quarks. The simple parametrisation:0 
x u v (x; q) oc x a (1 — x)@ , is not flexible enough to provide a good_pointwisc fit to 
our calculated distribution. However, the modernised fitting formliU 

xu v (x) oc x Vl (1 - x) m (1 - e u \fx + 7„x), (5) 

with parameters: r\\, r/2, e u , 7„, can describe our calculated result very well. 

The importance of this is apparent when one appreciates tlptjthe functional form 
xu v (x; q) oc x a (1 — x) 13 ', a = 0.67, (3 = 1.13, can be obtainedt3'tL3 via the evolution 
from qo — 0.35 GeV of u v (x) = 9(x)9(l — x), which latter distribution corresponds 
to the valence-quark carrying each and every fraction of the pion's momentum with 
equal probability, and to a momentum-independent pion Bethe-Salpeter amplitude; 
i.e., it is equivalent to representing the pion as a point particle. The convex- up 
character of our result is a characteristic feature of calculations in which the pion 
is described by a finite-size Bethe-Salpeter amplitude (see, e.g., Refs. Hence 
the convexity of the valence-quark distribution can be interpreted as a signature of 
binding in a nonpointlike pion. 
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While the moments of our calculated distribution are indistinguishable from 
those of that fitted to data, there is a pointwise discrepancy between our result 
and the data. Currently we judge that this discrepancy can be attributed to the 
restricted function space used thus far in parametrising pion data, especially as the 
modernised fitting form, Eq. (|5[), can describe our calculation. We would be much 
interested in a rcanalysis of existing data using the updated parametrisations and, 
indeed, in new experiments; e.g, Ref. fl3|| , with small errors on the valence-quark 
domain, which might expose what we have described as the signature of quark- 
antiquark binding in QCD's Goldstone mode. 
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